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Abstract 

We investigate the quantum integrability of the Landau-Lifshitz model and solve the long-standing 
problem of finding the local quantum Hamiltonian for the arbitrary n-particle sector. The particular 
difficulty of the LL model quantization, which arises due to the ill-defined operator product, is dealt 
with by simultaneously regularizing the operator product, and constructing the self-adjoint extensions 
of a very particular structure. The diagonalizibility difficulties of the Hamiltonian of the LL model, due 
to the highly singular nature of the quantum-mechanical Ifamiltonian, are also resolved in our method 
for the arbitrary n-particle sector. We explicitly demonstrate the consistency of our construction with 
the quantum inverse scattering method due to Sklyanin [1], and give a prescription to systematically 
construct the general solution, which explains and generalizes the puzzling results of [1] for the particular 
two-particle sector case. Moreover, we demonstrate the S-matrix factorization and show that it is a 
consequence of the discontinuity conditions on the functions involved in the construction of the self- 
adjoint extensions. 
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1 Introduction 

The Landau-Lifshitz (LL) model has been the subject of great interest in low-dimensional condensed matter 
physics as a model describing continuous classical magnets (for a review see [2,3]). In recent years, there 
was a surge of interest towards the LL model in relation with the gauge/string duality, where the LL model 
appeared on both sides of the correspondence [4-10]. In both one-dimensional magnetism in condensed 
matter physics, and in the context of the gauge/string duality it has become clear that the integrability 
plays the crucial role (for a review see [11-14]), allowing to construct the exact solutions and revealing the 
rich structure of the spectrum. Despite many years of investigation, only the classical theory of the LL model 
has been thoroughly investigated. The classical solitonic solutions have been found and discussed extensively 
in [15-20]. More complete classical analysis became possible after the classical integrability was established 
for the isotropic case first in [15,21], and for the general anisotropic case in [20,22]. The action-angle 
variables were constructed in [16], and in [23] the classical equivalence between the LL and the non-linear 
Schrodinger (NLS) models was established, relating the flat currents of the corresponding models by a gauge 
transformation. The quasi-classical spectrum was analyzed in [24] and subsequently in [17]. 

In contrast, the development of the quantum theory of the LL model was affected by a number of missed 
subtleties and nuances, which, as a consequence, led to the wrong quantization procedure and incorrect 
results [1,25]. Let us remind the general procedure of quantizing continuous integrable models. Apart from 
a few specific models, for example, the non-linear Schrodinger and the fermionic Thirring models, which 
can be quantized directly in the continuous case by means of the inverse scattering method as well as by 
coordinate Bethe ansatz, it is a standard procedure to consider first the lattice version of the continuous 
theory. This is done to deal with the ultraviolet divergences and regularize ill-defined operator product at 
the same point. The quantum Hamiltonian, and other conserved charges, can be found then by using the 
well-defined trace identities. Although there are many lattice models corresponding to the same continuous 
theory [26-31], the requirement of integrability restricts (but not eliminates completely) the choice of the 
corresponding lattice model. A systematic method that in principle is applicable for any integrable continuous 
model was outlined in [26] (see also [31]). In [26] this program has been successfully implemented for the 
NLS model and the sine-Gordon models. The difficulties with this procedure were already emphasized in the 
original paper. Even for simplest NLS model, the construction turned out to be quite non-trivial, and the 
resulting quantum Hamiltonian had a form describing interaction of eight nearest neighbors. For the sine- 
Gordon model, and in general for other continuous models, the quantum Hamiltonian is non-local, namely, it 
describes the interaction which depends on all lattice sites. An alternative method suggested in [27,28] states 
the existence of local quantum Hamiltonians for continuous integrable models, but its practical construction, 
based on the representations of the Sklyanin algebras [32,33], is a complex and in general unresolved problem 
(see also [34]). Although for the specific LL model, its lattice version is known, and can be obtained from 
the XYZ spin chain, for other more complex continuous integrable models, the resulting lattice regularized 
quantum Hamiltonian will in general be of a quite complex form, if constructed following the procedure 
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of [26] . The non-locality of the Hamiltonian may be a serious barrier to deal with, if one is interested in 
other subtle properties of the system. Thus, the quantization of the continuous integrable systems is highly 
desirable to carry out directly, without first passing to the lattice version. 

In [1] such program was first initiated, revealing a number of interesting questions and nuances, associated 
with quantization of the continuous integrable systems. There are several important points that should 
be addressed. To begin with, the usual method of constructing the quantum Hamiltonian and the other 
conserved charges, which works in lattice models due to the well-defined expressions and a few continuous 
models (the NLS model), does not work for the LL and in general for the majority of continuous models. 
The formal usage of the trace identities analogous to the ones of the NLS model or the lattice models, 
leads to wrong results [25]. There is, essentially, no effective or systematic method of constriicting the local 
conserved charges in continuous integrable models, as the quantum corrections modify the formal ill-defined 
expressions, which follow from the formal trace identities. In Sklyanin's original paper only the action of 
the two-particle quantum Hamiltonian was found, which involved some g\iess-work and consistency with the 
classical and quasi-classical cases. The constructions also required an unusual space of quantum states in the 
quantum-mechanical picture and very specific continuity properties of the functions involved. The higher 
n-particlc sector quantum Hamiltonian, its action and the quantum states have been unknown until now. 
Surprisingly, this matter has not been investigated in detail, despite its obvious importance. 

There are other interesting subtleties arising in the quantization of the LL model. Namely, in the 
anisotropic case the standard passing from the classical to the quantum transfer matrix does not work. 
Instead, the i?-matrix, as well as the monodromy matrix, are essentially guessed to satisfy the Yang-Baxter 
and bilinear relations. The construction of the monodromy matrix requires an additional spin operator, 
which, as a result, changes the algebra of the spin operators, thus, giving rise to the Sklyanin algebra [32,33]. 
Although the existence of such algebras was known before for the lattice systems, in the context of the 
continuous model its appearance is not very clear and remains an open problem. 

More importantly, it was realized in [1] that there arc essentially two distinct classes of the LL model, 
corresponding to the su{l, 1) hyperboloid and su{2) sphere cases. As it was correctly pointed out by Sklyanin, 
and missed by others who attempted to quantize the LL model, only in the su{l, 1) case one may construct 
physically meaningful states. In the su{2) case, the scalar product turns out to be not positively defined. 
Let us note, that the statement above is for the ferromagnetic case. It is possible to construct a positively 
defined scalar product for non-ferromagnentic vacuum in the su{2) case [35]. It is well-known, however, that 
in the case when the ferromagnetic vacuum does not exist, the algebraic Bethe ansatz is not applicable, and 
a more sophisticated construction, similar to the construction for the sinh-Gordon model, is needed [36,37]. 
This problem will be considered separately. 

In this paper we consider the isotropic su{l, 1) LL model and construct the desired quantum Hamiltonian 
for arbitrary n-particle sector, and its action on the states, which we also construct in detail. We recover the 
correct spectrum, and consider in detail the continuity properties of the functions involved. For the specific 
case of the two-particle sector, our results exactly reproduce Sklyanin's construction. The main point of our 
method is the regularization of the Hamiltonian directly in the continuous case. We achieve this by employing 
the split-point regularization, which effectively makes the Hamiltonian non-local, and constructing the space 
of quantum states, which requires a careful analysis of self-adjoint extensions in agreement with the scalar 
product. We emphasize that after removing the regularization, the Hamiltonian is local. As we show below, 
the quantum mechanical Hamiltonian for the LL model yields a highly singular potential which contains 
second derivatives of the delta-function. It is well known that in the space of functions with the usual scalar 
product it is possible to construct the self-adjoint extensions only up to the first derivative of the delta- 
function [38] . More precisely, there exists a 4-parameter solution in the space of function, that contains the 
derivative of the delta-function potential, and no solution exists which contains the second derivative of the 
delta-function. Thus, to circumvent this barrier, we propose to construct the self-adjoint extensions in the 
space of vector-like states, which is a novel feature, not considered in the literature previously, to the best 
of our knowledge. This provides an explanation to the ad hoc solution found by Sklyanin for the qiiantum 
mechanical Hamiltonian for the two-particle sector. Although we consider here in detail only the LL model 
as the simplest illustrative example, the above-mentioned difficulties and subtleties of a similar to the LL 
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model character will exhibit themselves also in other continuous integrable models. Thus, our method is 
general enough to be applicable to a wide-range of continuous integrable models. 

We also demonstrate in our method the S-matrix factorization, which is the underlying property of 
quantum integrability. It is worth noting here, that we have already considered the three-particle S-matrix 
factorization for the LL model in our previous paper [39] . There, we had shown that in the first non-trivial 
order of perturbation series, the three-particle S-matrix is indeed factorizable into a product of the two- 
particle S-matrices. Our work was based on earlier calculations of [40], where the LL model was considered 
from the field theoretic point of view, and the exact two-particle S-matrix was found by summing up the 
bubble diagrams, surviving in the two-particle scattering process (see also [41] regarding the problems with 
diagonahzation). We would hke to emphasize the difference between the su(l, 1) LL model considered here 
and the LL model considered in [40]. Only in the su(l,l) case one can construct a consistent quantum 
theory with the ferromagnetic vacuum. It appears that the LL model considered in [40] actually corresponds 
to the su{2) case. As we explained before, the construction of positive defined metric in the space of states 
corresponding to such a vacuum is mathematically challenging task, and was considered in [35]. The types 
of excitations in the su(l, 1) LL model and the case of [40] are also different. We do not see direct connection 
between the two models at this point, although this is an interesting problem to reconstruct the results and 
the spectrum of the [40] following the inverse scattering method and the method we propose in this paper of 
regularization of the local conserved charges together with construction of the self-adjoint extensions. It is 
not surprising then, that the n-particle S-matrix we find here for the su(l, 1) LL model (which is in complete 
agreement with the Sklyanin's result) is different, albeit by a coefficient, from the S-matrix found by [40]. 

Our paper is organized as follows. In section 2, we brieffy review the Landau-Lifshitz model in the context 
of the inverse scattering method. In section 3.1, we illustrate our method on the simplest two-particle case, 
give the regularized continuous quantum Hamiltonian, construct the self-adjoint extensions, and derive, in 
complete agreement with [1], the spectrum and the continuity properties of the functions involved. In section 
3.2, we consider the general n-particle case and show that the regularized quantum Hamiltonian used in the 
two-particle case is enough to construct the self-adjoint extensions and the spectrum in this general case. In 
section 4, we show the S-matrix factorization as the consequence of our construction. In section 5, we give 
a brief summary of our results and outline future problems. 



2 The LL model: Quantum inverse scattering method 

In this section we review, following [1] , the main features of the Landau-Lifshitz model in the context of the 
inverse scattering method, and discuss the arising difficulties and subtleties of this approach to quantization 
of the system. 

The Hamiltonian for the anisotropic Landau-Lifshitz model has the following form: 



(a,S,9,S) + 472 (5^)^-1 (1) 



where the vector S = {S\S'^,S^), and the scalar product is defined as (S, S) = (S^) - e (S^) - e {S^ = 1- 
The Poisson structure has the form: 

{SHx),S^{y)} = ±iS^{x)5{x-y) 

(2) 

{S-{x),S+{y)}^2ieS'{x)6{x-y) 

where 5*^ = S^±iS^. Here 7 is the anisotropy parameter, and the choice e = ±1 corresponds to the sm(1, 1) 
and su{2) cases correspondingly^ As explained in [1] the isotropic (7 = 0) and anisotropic (7 7^ 0) cases 



^As we mentioned in introduction, in this article we will only consider the su{l, 1) case, which, unlike the su{2) case, 
corresponds to the physically meaningful states in the ferromagnetic vacuum, with the particular choice of the representation 
for the states. 



4 



are essentially different, and should be considered separately. In the former case, the conventional inverse 
scattering procedure goes through without any changes - the Yang-Baxter and bilinear equations are satisfied 
with the appropriate choice of the i?-matrix. This is in contrast to the latter case, where the monodromy 
matrix and the spin operator algebra (Sklyanin algebra) have to be modified by hand for the intertwining 
relation (O to have a solution [32,33], with the i?-matrix corresponding to the XX Z model. The Sklyanin 
algebra naturally appears in lattice systems [32] as a consistency condition of the intertwining relations ^ . 
It is less clear, however, how to derive the Sklyanin algebra directly in the continuous case, and, in particular, 
in the LL model. We will consider here only the isotropic case. 

The Poisson structure ^ is replaced by the commutation relations for the S'-operators in the standard 
manner: 

[S^{xlS^{y)\ =±S^{x)5{x-y) 

(3) 

[S-{x),S+{v)\ ^2eS\x)5{x-y) 
The vacuum considered here corresponds to the ferromagnetic case: 

S\xm^e\Q) 

(4) 

S-{x)\Q) = 

The quantum £-operator in the isotropic case takes the form: 

C{u,x) = -[ )' 3)/ (5) 



u 



S-{x) ^S^ix) 



and the corresponding monodromy matrix is given by the expression: 

Here u is the spectral parameter and P signifies the path-ordered exponential. In this case the bilinear 
relatioiH 

(1) (2) (2) (1) 

R{ui - U2) T{ui)T {U2) = T {U2) T [ui)R{ui - U2) (7) 
is satisfied with the quantum i?-matrix given by the following form: 

R{u) ^ Wa{u)aa® (Ja (8) 

where the summation over the index a = 0, 1, 2, 3; wq{u) = u — i/2; wi{u) = W2{u) = ws{u) = —i/2. 

To construct the representations of ([3]) in the ferromagnetic vacuum, one writes the vector in the form, 
analogous to other continuous integrable models [31]: 

|/„) = / dxi...dXnfn{xi..Xn)S^{xi)...S+{Xn)\0) (9) 



where fii{xi x^) are continuous and decreasing sufficiently fast functions for the integral © to be well 
defined. A simple calculation shows [1] (see also the Eq. (HT]) below) that the scalar product {gn\fn) 
is positively defined only for the su(l, 1) case, while for the su{2) case, the matrix element is indefinite. 
Thus, only in the su(l,l) case one is able to construct physically meaningful states in the ferromagnetic 



(1) (2) 
''We use the standard notation T = T ® 1 and T 
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vacuum. Therefore, one has to take e = 1 in ([TJJ]). After passing to the infinite interval, the operators in 
the monodromy matrix ^ wiU satisfy the standard commutation relations, and together with the choice of 
the ferromagnetic vacuum, one can apply the well-known procedure of the algebraic Bethe ansatz to derive 
the spectrum and the eigenfunctions, which have the form (see [1] for complete details): 

\ui...un) ^ B{ui)...B{un)\0) (10) 

In the classical and lattice models, the bilinear relation ([7]) guarantees the existence of the integrals of motion, 
which are obtained from the generating functional I(u) = Tr [r(M)] . Here, however, one faces a difficulty, 
which is not present in the classical counterpart or the lattice version. Namely, in the classical case, one can 
simply decompose I{u) in the series : 

I{u)=J2hu'^ (11) 

k 

to obtain the local integral of motion. For example, the classical Hamiltonian ([T]) can be shown to be one 
of the charges in the Ik series. In the quantum case, however, the difficulty is that the local charges, in 
particular the Hamiltonian ([T]), contain operator product at the same point, thus making the local integrals 
of motion not well-defined quantities. Thus, the formally defined series ([TT|) cannot be used to obtain the 
integrals of motion, and, in particular, the quantum Hamiltonian is not a priori known. Thus, construction 
of the local integrals of motion turns out to be a non-trivial problem in the continuous quantum theory. 
From the field theory point of view, this should correspond to the renormalization procedure, which to the 
best of our knowledge has not yet been performed for the Landau-Lifshitz model. 

In [1] the quantum Hamiltonian was not found, and only the action of the (local) quantum-mechanical 
Hamiltonian on one- and two-particle sectors was presented. To formulate it, it was necessary to introduce 
new bosonic fields \l/„(a;), corresponding to the n-particle clusters, so that \l/„(a;)|0) = 0, satisfying the 
following algebra: 

[^Ux),^+iy)]^^nuA^-y) (12) 

One can then represent the S-operators as the following ?^-particle cluster decomposition: 

oo 

oo 

S+ix) = S+^+ix) + J2sl^t+li^)'fn{x) (13) 

oo 

S-{x) = so^iix) + ^s+'^'+(x)*„+i(x) 

71=1 



where Sq = 1; S([ = v^; = n; s+ = \/in+ l)n {n > 1). Using this cluster decomposition one can show 
that the two-particle eigenstate (|TOl) has the form: 



|U1,U2) - 



|0) 



2 cos Ml cos U2 

+ [ dxidx2 (c(pi,p2)e'(P^"^^+P^^^^ +c{pi^p2)e'^P''''+P''''^^ ^'+(xi)*+(a;2) 



where 



C{P1,P2) = TTf N 15 

2(pi -P2) 

Note, that in other known continuous models, solved by the coordinate Bethe ansatz, the first term in (|14p is 
absent. This is the case, for example, for the bosonic non-linear Schrodinger and fermionic massive Thirring 
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models. This somewhat unusual feature of the LL model, which was also discussed in [39], will be given 
explanation in the next section when constructing the self-adjoint extensions. 

As we mentioned earlier, the quantum field-theoretic Hamiltonian was not found in [1], but the action 
of the quantum-mechanical Hamiltonian on the two-particle sector was essentially guessed. To write it 
explicitly, it was necessary to introduce the space spanned by the states of the following type: 



l/> 



fi{x) 



(16) 



dxidX2f2{xi,X2)^t{xi)'^^{x2) 



|0) 



so that the function f2{xi,X2) is symmetric and smooth everywhere except on xi ^ X2 line, and fi{x) = 
f2{x,x). The Hamiltonian action is defined as follows: 



2{d,,~d,,)f2\:iz:iti-d^j2{x,x) 

- (^xi +di,) h{xi,X2) 



(17) 



It is not difficult to check that the above action of the Hamiltonian H2 is Hermitian with respect to the 
following scalar product: 

(18) 



\f\\'^ = J J dx\fif + J dxidx2\f2{xi,X2) 



Xx>X2 

One can check, that the solution to (|16m8p leads exactly to the two-particle state (|14l) . 

Although this form of the action was guessed by Sklyanin for the two-particle sector, its origin and the 
general n-particle Hamiltonian action and the corresponding space of states were unknown. In the next 
section both problems will be resolved, the local quantum Hamiltonian will be proposed, and the quantum 
states will be constructed. 



3 Quantum Hamiltonian and self-adjoint extensions 

As we discussed in the introduction, the construction of the local quantum Hamiltonian is a complicated 
task in general. The standard procedure of putting the continuous theory on the lattice to regularize the 
ultraviolet divergences leads generally to non-local Hamiltonians, and only the existence of the local form 
can be proven, while construction in practice is a complicated and unresolved problem. On the other hand, 
only a few continuous integrable models (NLS, massive Thirring model) allow direct quantization of the 
system by coordinate Bethe ansatz, without using the inverse scattering method. The essential difference of 
the LL model from the NLS or Thirring models is the presence of more severe singularities in the quantum 
mechanical Hamiltonian. Indeed, in the NLS model, the quantum mechanical interaction is described by 
the d{x) potential, while in the case of the LL model the interaction is highly singular and is proportional 
to dxdyd{x — y). The standard procedure to deal correctly with such singular potentials is to construct 
self-adjoint extensions. Even though for the NLS and Thirring models the problem was solved without 
constructing self-adjoint extensions, in general this is not correct. 

The construction of self-adjoint extensions for the LL model is, however, immediately bounded by the 
following fact. It is known that in the space of functions with the usual scalar product it is possible to 
construct the self-adjoint extensions only up to the first derivative of the delta- function [38], namely, there 
exists a four-parameter extension in the space of function, that contains the derivative of the delta-function 
potential. Since for the LL model the interaction is of the second order derivative of the delta-function, the 
above statement means that one has to construct a different scalar product in the new space of function. In 
the next section we will present this construction in detail first for the two-particle case, before generalizing 
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our analysis for the n-particle case. We will derive the Sklyanin's result pGlllSp and find its n-particle 
extensions. To do this, we will also propose the regularized continuous quantum Hamiltonian that will be 
checked to give correct results for any n-particle sector. 



3.1 Two-particle sector 

It is easier to demonstrate the idea of our method on the two-particle case, and in the more general case, 
considered in the next section, the complications are only of the technical character. Let us begin by writing 
down the regularized quantum Hamiltonian, corresponding to ([T]). It is not difficult to see, that the direct 
application of the Hamiltonian ^ to the two-particle state © 

1/2) - / dxdyh{x,y)S+{x)S+{y)\0) (19) 



leads to undefined singular expressions of the type d^6{x)\x=o- The idea is to regularize the continuous 
Hamiltonian by the split-point method. Namely, we take our quantum Hamiltonian to be of the form: 

Hq = limiJe (20) 



He^l I dudvF,{u, v) [-duS^dyS^ + duS+d^S- + d^dy {S^{u)S{u - v)) - dudj{u - v)] (21) 



where 

2 , 

Here the function F{u,v) is any smooth function, depending on some parameter e, so that 

limF,{u,v) = S{u~v) (22) 

We assume that F{u, v) decreases rapidly enough to make ((2T|) well-defined. Let us emphasize, that although 
the intermediate Hamiltonian (|2T|) is essentially non-local, we will remove the regularization ([22|) only after 
computations, thus, making the theory local. The first two terms in the limit ([22|1 will go to the classical 
Hamiltonian ([1]), and the last term of (PT|) is introduced to remove the infinities □ It is not difficult to show 
the following properties of the Hamiltonian Hq: 

Hq\0)=0 (23) 

Hglh) = J dx[-dl]h{x)S+{xm 

where is the one-particle state 

1/1) = / dxf,{x)S+{x)\0) (24) 



Thus, for the one-particle state we recover the correct solution fi{x) ^ e'x.p{ipx), with the energy Ei — p^. 
Let us consider now the more complex two-particle case. One can show that 



HQ\f2) = - y dxdy{dl+dl) h{x,y)S+{x)S+{y)\0) (25) 

+ J dx[id,-dy)Mx,y)\:zl+l-d.dyf2ix,y)U^y] S+ix)S+{x)\0) 

Thus, for p9p to be an eigenstate of Hq, we must require the following matching condition: 

{dx - dy) f2{x,y)\lzl~l - dxdyf2{x,y%=y = (26) 



^Note, that this term is enough to remove infinities in all n-particle sectors simultaneously (see the Eq. I I27I I). 
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We note here, that in the more general case for the n-particle sector, these matching conditions, resulting 
after the action of the quantum Hamiltonian (|20p on the n-particle state, will remain the same, with the 
obvious change x ^ Xi and y ^ Xj. After some straightforward algebra, we can derive the action of the 
quantum Hamiltonian (|20p on the n-particle state: 



Hq\U) = - / dl?(A/(l?))[]5+(x,)|0) (27) 

i—1 

From here we immediately derive the general matching conditions: 

- {d,f{x) - ddml';Zt\ = dAf{^%.=-, > > J (28) 

This is, in fact, the reason for the S-matrix factorization, as we will see below. We also emphasize, that the 
quantum Hamiltonian (|20p does not acquire any further corrections in the higher n-particle sectors. 
We will show now, that the matching condition (|26p together with the equation following from (pS)) 

-{dl+dl)f2{x,y) = {E2)f2{x,y) (29) 



where E2 is the energy of the two-particle state, recovers the Sklyanin's solution (|T4|). We will first construct 
the space on which the quantum-mechanical Hamiltonian acts. 
Let us consider a space V generated by the vectors of the form 

* - i Itl ) (30) 



where the function fi{x) is determined by f2{x,y) and, possibly, its derivatives &i x ^ y. The actual form 
of fi{x) will be fixed later. For a given non-negative number, a, we define a scalar product on V as follows: 

00 

j gl{x)h{x) dx + JJ g;{x,y)f2{x,y) dxdy (31) 
—00 

where 

/f\ ) and$=f ) (32) 

V Mx,y) J V 52(2;,?;) J 

We require the function fi{x) to belong to £^(R, dx) and f2{x) to £^(M^/{x ~ y}, dxdy). Further conditions 
on fi{x) and f2{x,y) will be imposed later. 

Now we define the operator, H, on V with the following properties: 

i) acting on f2{x) it is simply the Laplacian —A = —d^ — dy everywhere in R^/{x = y}, i.e. 

H^=( , ) (33) 



-A/2(a;,y) 

with some operator h to be determined later; 

ii) it is Hermitian with respect to the scalar product ([?T|) . i.e. 

= (34) 

Using 

9*{x,y)l\f{x,y) = ^^{g* {x,y)d^f{x,y)) - ^^{{^^g{x,y)y f{x,y)) + {Ag{x,y)y f{x,y) 
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where i = {x, y}, and not assuming continuity neither functions nor their derivatives at a; = ?/, we have 



92ix,y)^f2ix,y) dxdy^ JJ {Ag2ix,y)y f2{x,y) dxdy 
dy [9*2i.x,y)dxhi.x,y) - {dxg2{x,y))* f2{x,y)]lzltl 



dx [g2{x,y)dyf2{x,y) - {dyg2{x,y))*f2{x,y)]lJltl (35) 

Now, to be more specific, we are going to impose some conditions on f2{x,y). We require it to be 
continuous at a; = ?/0 In this case Eq. (j35|) simplifies to become 

92{x,y)/^f2{x,y) dxdy ^ j j {Ag2{x,y))* f2{x,y) dxdy 

x^y x^y 

oo 

dx 152(2;, a:^) [dxf2{x,y) - dyf2{x,y)fyllt\ 

-[{d^g2{x,y)r - (9,52(2;, y))i^::+:/2(x,x)} (36) 

ft is obvious that the first line in (|36p has needed for Hermiticity form while the second and third lines 
should be compensated by h (see Eq. (p3| ) in such a way that the full H would become Hermitian. Before 
proceeding, we must fix the relation between fi{x) and f2{x,y). Looking at the one-dimensional integral 
in p4|) . one sees that a natural choice is fi{x) = f2{x,x) (which is possible after we required continuity of 
f2{x,y)). Now it is not difficult to see that the following form of h 

hhix) := i [d.j2{x,y) - dyf2{x,y)]lZltl + hh{x) , (37) 

where hi is any Hermitian in £^(R, da;), does the job. 

Thus, we see that the requirement for H to be Hermitian (or rather symmetric) with respect to the 
scalar product (|3ip does not fix the Hamiltonian completely even after we had chosen some conditions on 
the components of |^). But, in fact, we still have to use one consistency condition: the image of I^P) under 
the action of H should belong to the same class of vectors, namely the first component should be related to 
the second one 

hh{x) := - [dxf2{x,y)^dyf2{x,y)]lZltl + hihix) = ^Af2{x,y)U=y (38) 

This put some constraints on the form of hi as well as imposes some conditions on the behavior of f2{x, y) 
a.t X = y. 

For the LL model, it is not difficult to show from (fTO|) . that the coefficient in the scalar product PT|) 
a = 1/2, and using the Eq. (|26p we find that hi = —d^. Thus, we obtain from ([55)1 : 

/ 2id,,-d,,)f2\:iz:\ti-d'j2{x,x) 

H2\.f) = I (39) 

V -{dl,+di;)f2{xi,X2) 

This is exactly the formula pT|) guessed by Sklyanin. 



*The other possibilities are also interesting but we consider the one that is relevant for our problem. 
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3.2 n-particle sector 

Let us consider a space of vectors of the form 



|/„) = / d"f fix) S+{xi) ■ ■ ■ S+{Xn)\0) (40) 



Then it is not difhcuh to calculate a 'scalar product' in this space. Let {X^} be a partition of the set {xi}: 

Mp 

U ^™ = {^0 ^^'^ ^™ Pi -'^n = SrnnXm 

and let tm be a 'collective' coordinate for all Xi G Xm- Then the resulting 'scalar product' is 

(ff«|/n>= E ^''-"'Cp j d''-tg*{x)f{x)\{,^^x„.^^u. (41) 

partitions 

Here Cp are some combinatorial factors, which are positive. From here it is obvious that this 'scalar product' 
will be an actual scalar product only for the non-compact case of suil, 1), where e = 1. From now on we 
concentrate on this case postponing consideration of the su{2) model to future work. It is instructive to 
write down n — 2 and n = 3 cases explicitly: 



n = 2: (.g2|/2>-4 
" = 3: (.g3|./3)=8 



2 j dxdy g*{x,y)f{x,y) + J dx g*{x,x)f{x,x) 
6 / dxdydz g* {x,y, z)f{x,y, z) 



(42) 



+ (3-3) J dxdy g*{x,x,y)f{x,x,y) + J dx g*{x,x,x)f{x,x,x) (43) 

While (dH) is exactly the scalar product used in [1] (see also the Eq. ([3T|) ). the formula is presented 
here to explicitly demonstrate the first non-trivial case. 

There are a couple of useful interpretations of (|4T|) . The first one is along the lines of [1]. Namely, one 
can think of (f40|) as a vector-function with (f4T|) as a natural scalar product. Even though this is a useful 
interpretation, mostly because this nicely fits into the cluster picture of [1], the more mathematically rigorous 
one is as follows. Namely, it is a scalar product in L^(i?), where the completeness is defined with the help of 
the Riemann-Stieltjes integral, rather than just the Riemann one. Here, /i is the 'measure' for this integral, 
which formally could be written as 

' f{x)d"fi{x) , (44) 

where fj, now is not required to be a smooth function. The only requirement is that / and /i do not have 
discontinuity at the same points, which is the case - our function is continuous. And that is why it is that 
hard to define the self-adjoint operator - one needs to be extremely careful at the points, where the measure 
is discontinuous. The measure fi is very easy to read from the scalar product. In particular, for the case 
n = 2, Eq. ((42|) . we have (formally): 

y) = i/2xy - x6{x -y)- yO{y - x) (45) 

Regardless of the interpretation, the strategy in defining a self-adjoint operator is the same: because 
the continuity is a property only of the function itself but not of its derivatives, there will be surface terms 
that will compete with lower dimensional integrals. They should be accurately taken into account. That is 
exactly what we do below (and what we saw above for the n = 2 case). 
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Then we proceed exactly as we did in n = 2 case: we will construct a self-adjoint extension of the 
Hamiltonian, which in the 'bulk', i.e. everywhere except Xi = Xj for all possible i and j, reduces to a 
Laplacian: 



N 



E^' (46) 



i=l 



As usual, this will amount to imposing some sewing conditions on derivatives of f{x) at Xi = Xj. Let us 
start with the integration of the highest dimensionality in (|4ip . i.e. when all the clusters contain just one 
coordinate, Xi = Xi. Performing integration by parts and taking into account surface terms, we have 



d''xg*{x)dff{x) = J d^x {dfgnx)f{x)+ 

Summing over i, we have the result for the full Laplacian. It is clear that the action of the Hamiltonian on 
{n — l)-components of the aforementioned vector-function, i.e. when only two of the coordinates are equal, 
should exactly compensate the unwanted term 

EE/ II'^-'^'^ (48) 

and simultaneously reproduce the term 

- E E / n [3*(-)5./(x)]::::^:;: (49) 

What is the Hamiltonian acting on fij :— f{x)\xi=xj that does this? The answer is obvious, if we recall that 
the function / itself is continuous. Then we can define 

kjh, = -« (aj(f) - d,f{x))l%l'+\ + . (50) 

Here a is a combinatorial coefficient expressed in terms of Cp from (|¥T|l (in fact, it is not hard to see that it 
is always the same as for n — 2: a — 2, e.g., it is explicitly seen for n — i from (|43p : a = | = 2); hij is some 
operator in [n — l)-dimensional space, such that its non-Hermiticity is on the next, (n — 2)-dimensional, level. 
In a moment we will see that it is nothing but a Laplacian acting on fij . It is also easy to understand why 
there are two terms in (jSOp instead of the one, as one would nai'vly expect from (j49|) : the term at Xi — Xj 
will arise twice in the sum - the first time as 



[5*(x)aj(f)]^:z:^+: (si) 



and the second time as 



[9*imfim:2'r: (52) 



The residual freedom in the definition of hij coming from yet undefined hij cannot be fixed only by 
requiring Hermiticity - hij (and its analogs in lower dimensions) remain undetermined. It is removed by 
requiring that it should agree with the matching conditions that follow from the explicit action of the 
quantum Hamiltonian on the state (j40p . As we had already mentioned in the previous section, the matching 
conditions, following from the action of the quantum Hamiltonian (|20p on the n-particle state (|40p . are, 
essentially, the same as in the two-particle sector case (pS)) . Below we show that providing the matching 
conditions (|28|) 

- id, fix) - 5./(f)):^::;;: = dAf{^)U.=^, , Vz > J (53) 
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then hijfij is, in fact, equal to Af{x)\xi=xj if hij equals to (n — l)-dimensional Laplacian. 

Let us find how a Laplacian in a lower dimensional space is related to the one in n-dimensional evaluated 
when some coordinates coincide. For the future use we consider more general case when there is one cluster 
but not necessary of the size two, {xki, As above, let i be a collective coordinate for this cluster 

and i runs over the rest of the coordinates. Then the lower dimensional Laplacian is 

Afk,...k,{t,x.) = id^+J2df)fit,^^) (54) 

i 

Now for dff{t^ Xi) we have 

d^fit,x.)^l J2 ^^/(^) +2 E d,d,fix)\[x,}=t (55) 

and finally 

M{^)\{x^}=t^Mk^...h,[i,Xi)-2 ^ d,dj![x)\{x^^^t (56) 

Now we see that if we define hij as 

= -2 (a./(a-) - d,j{x)f:^^^ + A,.,/,, (57) 
and express Ai=jfij using ([56]), we will have 

k,f,,=Af{x)U=j (58) 

which, of course, should be the case as Hamiltonian should take continuous function to continuous function. 

We can repeat each step starting at any level (n — k) and using collective coordinates, ti. Because the 
non-Hermitian part on this level is in corresponding h, which has exactly the form of (n — l)-dimensional 
Laplacian, we will arrive at the level {n — k — I) with exactly the same form of the lower dimensional 
Hamiltonian. This completes the prove that the Hamiltonian, given as a Laplacian in the 'bulk' plus the 
matching condition (j53p . is, in fact, a self-adjoint operator in the space L^(i?). 



4 S-matrix factorization 

Although the S-matrix factorization is the underlying property of quantum integrable systems, it is quite 
difficult to analytically prove it using the standard perturbative calculations. Until now such calculations 
were fully performed in all orders only for the non- linear Schrodinger model [42,43], where the S-matrix 
factorization was indeed proven - first for the three-particle scattering process, and then generalized to 
the rt-particle scattering process. Recently, in [39] we have considered the S-matrix factorization following 
similar perturbative calculations, based on earlier two-particle S-matrix calculations of [40]. We were able 
to show the three-particle S-matrix factorization in the first non-trivial order in the perturbation series. 
As we explained in [39], there are essential differences between the S-matrix calculations for the NLS and 
LL models. Besides the technical difficulties associated with rapidly increasing, in each perturbation order, 
number of vertices in the LL model, which leads to a complex diagrammatic analysis, there are conceptual 
difficulties related to the identification within the field theoretic approach of the Bethe particles. 

Here, however, we will establish the S-matrix factorization of the su(l,l) LL model in a more direct 
fashion, following the results of the previous sections, where we have explicitly constructed the self-adjoint 
extensions and their continuity properties. Thus, having the exact expressions for the rt-particle wave- 
functions and using the matching conditions, it is not difficult to establish the S-matrix factorization. In 
fact, one can proceed in the same fashion as for the simpler NLS model (see for example [31]). Indeed, the 
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main result from the previous section, that aUows explicit calculation of the n-particle S-matrix, is that the 
matching conditions for the n-particle case (1551) are exactly the ones for the two-particle case (PB)) . In other 
words, it is enough to solve the equations for the two-particle case, in order to obtain the solution for the 
general n-particle case. This is, essentially, the S-matrix factorization. 

To construct explicit expressions, we follow the derivation of [31] for the NLS model. It is easy to see 
that for the LL model, the matching condition ([28|l leads to the following n-particle wave-function: 

fnix^\pJ) = (const)^(-l)We.?i *^*'']^ fp^.^ _ ^ip^^^p^A (59) 



{i} 



where {i} denotes all possible permutations of (1, ...,n). Thus, the n-particle wave-function is factorized in 
terms of the two-particle wave-functions. As a consequence, the n-particle scattering S-matrix has the form: 

Sn{pi, --.Pn)^ WS2{Pi,P]) (60) 

where the two-particle scattering S-matrix for the su{\, 1) LL model has the form: 

2{pi ~ P2) - ipiP2 



2(j3i - P2) + ipiP2 



(61) 



Let us note, that for the NLS model the expression for the wave-function has similar to ([55]) form, where 
the third term in the brackets is a constant. Here, we have momenta product instead, which is the result 
of derivatives present in the interaction vertex (see for details [40]). With this in mind, wc can intuitively 
think of the LL model as the NLS model with momentum-dependent interaction. 



5 Conclusion 

We have considered the quantum intcgrable properties of the Landau-Lifshitz model, and proposed a method 
to construct the quantum Hamiltonian. Most importantly, we achieve this directly in the continuous case 
by regularizing the ill-defined Hamiltonian, and constructing the necessary self-adjoint extensions. This 
method allowed us to consistently derive the spectrum, which we show to coincide with the one following 
from the quantum inverse scattering method. We gave an explanation and derived in the most general n- 
particle case the puzzling construction of Sklyanin [1] (for the particular two-particle sector) of the quantum- 
mechanical Hamiltonian action on the vector-like state. The continuity properties of the functions involved 
in the construction of such states have also been carefully investigated. These properties are defined by the 
corresponding matching conditions, and, as we have shown, lead to the S-matrix factorization property. 

The particular difficulties of the LL model quantization are the ill-defined operator product of the local 
conserved charges, as well as the highly singular potential in the quantum-mechanical picture, which make it 
impossible the use of the trace identities in the quantum case. Thus, it is clear, that the method considered 
in this paper should be applicable to any continuous integrable model which has a singular nature. Since 
we have only considered the isotropic LL model, the next natural problem is to consider the anisotropic 
LL model, which is of great importance in the theory of integrable models. This, however, seems to be a 
more complex task, as the algebra of observables should be, for consistency, modified by hand in the inverse 
scattering method, forming the Sklyanin algebra. Let us remind, that this algebra is naturally obtained from 
the lattice models, but its appearance in the continuous models is less clear. It would be interesting to give 
a direct derivation of the Sklyanin algebra without appealing to the lattice version, since in more complex 
continuous integrable models the construction of the corresponding lattice models, as we have discussed 
in the introduction, is generally a quite complex task, that has not been well-understood even for simple 
models. 

Another important problem, also discussed in introduction, is to investigate the su{2) LL model. Let us 
remind that in the Sklyanin's original paper [1] as well as in our work, only the quantization of the su(l, 1) 
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LL model is considered. This is done to have physically meaningful states in the chosen representation for 
the states, consistent with the ferromagnetic choice of the vacuum. Constructing representations in the su{2) 
case with positively defined metric seems to be a more complex task which has not so far been considered 
in connection with the quantization of the LL model. It is known that such representations are possible 
to construct, but the vacuum will not be ferromagnetic anymore. Thus, the algebraic Bethe ansatz is not 
applicable in this case, and one has to consider more sophisticated methods of finding algebraic solutions, 
much like for the sinh-Gordon model [36,37]. 

Although we have considered only regularization and diagonalization of the first non-trivial conserved 
charge (the Hamiltonian) , integrability implies the conservation of the infinite tower of charges, that should 
be possible to regularize in the manner similar to the method proposed in this paper. We do not currently 
know whether it is possible to do in a unified manner, or each charge should be considered separately. 
Clearly, the questions posed above for the LL mode, as the simplest representative of associated difficulties, 
will appear in other more interesting continuous integrable model. As an example, we mention the recently 
discovered Alday-Arutyunov-Frolov fermionic model, which appears in the su(l, 1) subsector of the AdS^ x 
strings [44]. There we expect the similar difficulties to appear in the quantization process, as the singular 
nature of the fermionic interaction terms will clearly require construction of the self-adjoint extensions 
and careful consideration of the conserved charges. These and other related problems are currently under 
investigation. 
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